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We discuss quantum dynamics in multi-dimensional non-linear systems. It is well-known that 
wave functions are localized in a kicked rotor model. However, coupling with other degrees of freedom 
breaks the localization. In order to clarify the difference, we describe the quantum dynamics by 
deterministic rigid trajectories, which are accompanied with the de Broglie-Bohm interpretation 
of quantum mechanics, instead of wave functions. A bundle of quantum trajectories are repulsive 
through quantum potential and flow never to go across each other. We show that, according to the 
degrees of freedom, this same property appears differently. 
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^ \ I. INTRODUCTION 



When the question "What is quantum chaos?" is thrown to us, we usually answer that no chaos is found in the 
quantum world. Successively, we cite a kicked rotor model |^ as a typical example. Its classical dynamics is reduced 
to the standard map, which is a well-known chaotic system. For a classical ensemble of the initial state, the dispersion 
psj ' of the momenta eternally continues to increase as a result of a random walk. For a quantized kicked rotor, however, 
^ " numerical simulation reveals that the diffusion in the momentum space begins to be saturated after a time scale, 
and the wave function is found to be finally localized Non-existence of quantum chaos can be also explained 

[ in general terms. For an isolated bound system, the energy spectrum is discrete and almost periodicity appears in 
. the evolution of the wave function even if the energy spectrum is coarse-grained . Further, even if the potential 
I ' is allowed to be temporally periodic, the wave function of the bound states also evolves almost periodically as long 
as the system is not resonant In the proofs of these theorems, the property that the Schrodinger equation is a 
^\ kind of wave equation and a wave function is a superposition of the oscillating solutions works as fundamental one to 
quantum mechanics. Therefore, quantum chaos is occasionally called quantum chaology or quantum psuedochaos |^ 
^ ^ . As reflection of this situation, the energy level statistics and the complex structure of eigenfunctions in quantum 

I non-linear systems rather than the dynamics have been searched for remnants of classical chaos. 
^ A decade ago, Adachi, Toda, and Ikeda |^ analyzed a system of coupled kick rotors which are classically chaotic. 

d ' They found that localization of wave functions is broken in the model, while it appears when the coupling between 
^ , the two rotors is absent. Then, Ikeda, Adachi, and Toda |^ showed that a quantum system with helper degrees of 
^ freedom monotonously absorbs energy from the external field [0. Recently, Kubotani, Okamura, and Sakagami |l|] 
J> ] have proposed that the coupled kicked rotors also yield quantum noise with decaying correlation and break irreversibly 
the quantum coherence in contacted other degrees of freedom. These results show that redundant degrees of freedom 
may drastically alter the quantum dynamics. Further, if we are allowed to apply the Schrodinger equation to the 
. whole universe, which is isolated and composed of infinitely many degrees of freedom, the fundamental questions 
" " ' arises: Why some degrees of freedom behave classically and others quantum-mechanically? Why only the formers are 
able to behave chaotically or irreversibly? These problems cannot be answered straightforwardly jl^ . Even without 
taking the ideal and mathematical limit of infinite degrees of freedom a priori, therefore, the quantum dynamics still 
involves important and interesting phenomena, although chaos in the exact sense, which is believed to contribute to 
irreversibility in the classical world [l^ , may not be present in the quantum world. 

According to the de Broglie-Bohm interpretation of quantum mechanics [ p^ , we can introduce a bundle of rigid 
trajectories which are consistently equivalent to a wave function. At first sight, the rigidity may be inconsistent with 
the uncertainty principle in quantum mechanics. Indeed, the uncertainty principle makes a product of the expected 
values of canonically conjugated variables superior to the Planck constant h. In the Copenhagen interpretation, the 
relation is considered as inaccuracy of the simultaneous measurement of the conjugate variables. However, there is no 
mathematical inconsistency between the quantum trajectory picture and the time evolution of a wave function driven 
by the Schrodinger equation as will be seen in Section 2. The de Broglie-Bohm interpretation considers quantum 
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expectation values as statistically averaged ones and thus the apparent inconsistency is due to the confusion of the 
measurement theories on which we base the interpretation of the uncertainty relation. 

Although the recovery of classical diffusion is confirmed by the previous works 0- [Q, it does not necessarily mean 
that the classical dynamics is restored in quantum mechanics. The dynamics needs to be estimated more directly. The 
associated problems to be resolved are how nonlinearity work in quantum systems and what role redundant degrees 
of freedom play. In order to reveal them, in this paper, we propose to utilize quantum trajectories. For analysis of 
wave functions, indeed, the Wigner and Husimi functions have been used p^ . These projective representations of 
wave functions on the phase space can be compared with classical statistical distributions. However, since they do 
not follow the Liouville equation, the comparison estimates only the correspondence or difference in the snap shot 
between the classical and quantum states. For the analysis of the dynamics, therefore, the description of quantum 
mechanics by trajectories is favorable, since the trajectories can be easily compared with classical orbits. By the 
quantum trajectory picture, even the standard notions of chaos, such as the Lyapnov index and the KS entropy are 
also introduced naturally into the quantum world p5| ] . 

This paper is organized as follows. In the section 2, wc introduce a coupled kicked rotator model as a multi- 
dimensional non-linear system and discuss the quantum picture of trajectories. We calculate the quantum trajectory 
of the system in the section 3. In the section 4, we give the summary and discuss the results. 



II. MULTI-DIMENSIONAL NON-LINEAR MODEL 



The kicked rotor model was first introduced by Chirikov [y. Its property has been widely investigated both in 
classical and quantum mechanics, since it has the dynamics which is typical of nonlinear systems. The Hamiltonian 
is 

oo 

Hk{q,p;t)^-p^ + kcos{q)Y,5{t-nT) (0 < g < 2^), (2.1) 

n=l 

where q and p are an angle variable with period 27r and its canonically conjugated action variable, respectively, and 
k and T are constant. The first term is quadratic with respect to the action variable p and shows a nonharmonic 
oscillator. Alternatively we can also consider it as the kinetic term with respect to the momentum p. The second 
one shows the temporally periodic and instantaneous kick, whose strength and period are parameterized by k and T, 
respectively. The classical dynamics is reduced to a discretized map. The variables (qmPn) = {q{t = nT + Q),p{t = 
nT + 0)) at the discretized time t = nT + (n = 0, 1, 2, • • •) are transformed as 

qn+l = qn + Tpn, 

Pn+1 ^Pn + fcsin((7„+i). (2.2) 

Here -1-0 indicates the time just after the kick. The transformation rule ( |2.2| ) is called the standard map, whose 
dynamics is characterized by one parameter K = kT. When if < 1, the variables {qn,Pn) are trapped on the KAM 
torus. When K exceeds the critical value Kcr = 0.9716- • •, they diffuse into the chaotic sea on the phase space. 

We couple two kicked rotors whose kick strengths are fci and fca. The total Hamiltonian of the coupled rotors is 
given by 

H{qi,pi,q2,P2;t) Hk,{qi,pi;t) -f ii^^ (ga , P2 ; + H„it{qi,Pi,q2,P2;t). (2.3) 

This type of model has been used by Adachi et al , Kubotani et al ||ll| , and Sakagami et al . In this paper we 
choose the interaction Hamiltonian Hi„t as 



Hint{qi,pi,q2,P2;t) = Cpppip2, (2.4) 



where Cpp is a coupling constant. The interaction Hamiltonian Hint, Eq. (2.4) causes continuous interaction between 
the two degrees of freedom through the action variables pi and p2- Hereafter we call the coupled kick rotors as a 
pp-coupling model. 



We quantize the coupled kicked rotors (p. 31). In the Shrodinger picture, we write down the time evolution operator 



U{t) = Toexp[-- / dsH{s)], (2.5) 
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where Tq denotes the time ordered product, and H{t) i s th e quantized Hamiltonian operator which is derived by 
replacing the canonical variables qi, q2, Pi, and p2 in Eq.(^.3|) with the corresponding quantum operators qi, (72, Pi, 
and p2. In general, incommutability between the conjugate operators in the Hamiltonian complicates estimation of 
the time ordered product. For the pp-coupli ng m odel, we only have to take account of the ordering at the time when 
the kick is added. As a result, the operator (2.5) is reduced to a simple form 

U{t) = Ui{h){U2Ui{T)r, (2.6) 

where 

Uiit) = exp[-^(p2 +pI + 2cpppip2)t], (2.7) 

U2 = exp[--(fcicos(gi) + fc2Cos(g2))]. (2.8) 

Here to = t — uqT and Uq = [t/T], where [x] denotes the maximum integer which is not superior to x. Note that 
the operators Ui and U2 are constructed only from the the action variable operators pi and p2 or the angle variable 
operators qi and q2, respectively, and can be expressed as diagonal matrices by choosing an appropriate representation. 
Therefore, the operation of Eq. ( |2.6|) to a wave function is decomposed into the two procedures: the interchange of 
the representation of the wave function and the multiplication of complex factors. The simplicity guarantees low 
roundoff error in the numerical simulation of a wave function. 

According to the de Broglie-Bohm interpretation of quantum theory , we construct a bundle of deterministic 
rigid trajectories from a wave function in the pp-coupling model. First, we express the wave function with respect to 
the polar coordinate: 

$(gi,g2;t) = R{qi,q2;t)exp[^S{qi,q2;t)], (2.9) 

where the qi, q2 representation is chosen, and R and S are real- valued functions. Using this expression, we decompose 
the Schrodinger equation into two parts: 

5 „9 d „,,dS . ,dS.. d ^r,,,dS . ,dS.. „ „x 

and 

m + 2 + 2 ^ "^^'^^''-^'^ "^^-^^'^^^^^ ^ ^'-''^ 

where V{qi^ q2] i) is a original potential: 

V{qi,q2;t) = (fcicos((Zi) + fc2Cos((Z2))5Z'^(^-"T), (2.12) 

n=l 

and Vq((3'i, 92) is what we call a quantum potential: 

^ ,d^R d^R „ d'^R ,,2 rr. ox 

V,(qi,q2;t) = --(^ + ^ + 2C,,^);.^ (2.13) 

Next, we consider an ensemble consisting of quantum particles with the initial distribution function R{qi, q2',t — 0)^. 
For the particles, we can define rigid trajectories as follows. The momentum {pi,p2) of a particle on a trajectory is 
equalized to the spatial gradient of the phase part of the wave function, S: 

P^ = i^) (^ = 1,2). (2.14) 

oqi 

From the Hamilton's canonical equation, the particle moves according to the velocity: 

qi=pi+ CppP2 = (-r— ) 4- Cpp( — ), 
oqi dq2 

dS dS 

q2^P2+CppPi^{ — ) + Cpp{ — ). (2.15) 
oq2 oqi 
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I nteg ration of Eq. ( 2.15| ) with respect to t yields a time-parameterized trajectory. The relation (2.15) allows Eq 



( ^.10 ) to be rewritten as 



+ -^qiR^ + -^hR" = 0, (2.16) 
at dqi dq2 



which reserves conservation of the distribution function R{qi, q2',t) ■ We also note that Eq. ( 2.11 ) corresponds to the 



Ham ilton- Jacobi equation in classical mechanics. Compared with the classical one, the additional term Vq appears in 



Eq. (2.11). The term shows quantum effect on quantum trajectories. That is why Vq is called a quantum potential. 



We mention that Eq. (2.14) appears to show WKB approximation. By definition, however, S is not equivalent to 
the phase part of the WKB wave function, S'wkb- •S'wkb is determined by the classical Hamilton- Jacobi equation 
which does not include the quantum potential Vq. Therefore, only if Vq is negligible, S happens to be equal to 
S'wkb- In other words, the quantum trajectory picture is not limited to WKB regions and always defined consistently 
irrespective of any approximation schemes. 

III. QUANTUM TRAJECTORY IN NON-LINEAR SYSTEM 

In this section, we numerically trace quantum trajectories in the coupled kick rotor model constructed in the section 
2 and compare them with classical o nes. For this aim, we first solve the original Schrodinger equation instead of the 



quantum Hamilton- Jacobi equation (2.11). To represent a wave function numerically, we use 4096 x 512 mesh points. 
Setting h = 2Tr X 43/4096 = 0.0659 • • •, we can describe the momentum eigenstates with the eigenvalues pi from — 437r 
to 437r and p2 from — 437r/8 to 437r/8. For the model to analyze, we fix the parameters as ki = 2.0, k2 = 0.9, and 
T = 1.0. As an initial wave function, we choose a product of the momentum eigenstates for each degree of freedom: 

$0(<Zl,g2) = $Pio(9l)® *P2o(92) 

\ 2 TT 1 Z TT 

exp(--pio(gi - t)) X -==exp(--p2o(92 - 7)), (3.1) 



where piQ and p2o are initial eigenvalues. In the following concrete calculation, we set pig — P2q — 7r/2 . The wave 
function at the arbitrary time is obtained by multiplying the initial wave function by the o perator (|^). From the 



time-dependent wave function, next, we estimate the velocity of a quantum particle, (2.15). The spatial derivative 
of the wave function just on the mesh point is estimated by the inverse Fourier transformation. The value between 
the mesh points is interpolated p^ . A quantum trajectory is calculated by integrating the velocity field by the 
Rungge-Kutta method. 

We begin by observing gross time evolution of quantum state with respect to the dispersion of the momentum 
pi. Fig. 1 shows the second order moment Q =< {pi— < pi >)^ >, where <> denotes the expectation value 
estimated with respect to the wave function at the time t — nT + 0. The lines (No) and (P-P) indicate the single 
kick rotor model (fci = 2.0, Cpp = 0) and the pp-coupling model (fci = 2.0, k2 = 0.9, and Cpp = 0.2), respectively. 
For the single rotor, localization of a wave function is found. For the coupled kicked rotors, on the other hand, the 
diffusion in the momentum space seems to continue. In Fig. 2, the distribution of the momentum pi at t = lOOT -I- 0, 
I < pi\^{t = lOOT-l-0) > p is shown for each model. The distribution (No) is approximate to ln| < pi|<l) > p cx — |pi|, 
although the classical law for the chaotic rotor would be ln/(p) esc —p^, where f{p) is a classical distribution function. 
Comparison with the quantum kicked rotor modulated by classical noise is discussed in |^ . 

Fig. 1 shows that introduction of coupling with other degrees of freedom may change the quantum dynamics 
drastically. However, the recovery of classical diffusion found here does not necessarily mean that the Liouville 
dynamics is restored in the evolution of the wave function. Fig. 2 shows the plateau with 2it period appears in 
the distribution for the coupled rotors. As well as Fig. 1, it suggests similarity between the classical and quantum 
dynamics 

The classical dynamics of a single kicked rotor is reduced to a map between the variables at discretized time as 



mentioned in the section 2. For the coupled kicked rotors (2.3), the map is 



qi,n+l = qi,n + TPi.n + CppPj^n (i = 1, 2, j ^ i) , 

Pi.n+1 = Pi,n + fciSin(qi,„+i), (3.2) 

where {qi^n,Pi,n) — {qi{t — nT + 0),pi(t = nT -|- 0)) are the canonical variables at the discretized time t = nT + 
(n=0, 1, 2, • • •). The classical particle in the region [0, tt] is accelerated and one in the region [7r,27r] is decelerated by 
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a kick. Until the next kick, the momentum is conserved and the velocity of the particle is constant. The dynamics 
results in effective stretching and folding of a phase space volume element. 

For the quantum particle, on the other hand, the velocity can be altered, even if the original potential in the 
Hamiltonian is not present. The instantaneous kick yields the inhomogeneous velocity field of the quantum particle. 
After the kick, as a result, the particles are gathered inhomogeneously and high density regio ns are produced. The 



gradient of the density of the particles gives nonzero quantum potential through Eq. ( 2.13 ). In other words, the 



quantum partic le feels repulsive force from the other particles. Further, uniqueness of the solution of the differential 



equation (2.15) guarantees that two quantum trajectories cannot go across each other. Fig. 3 shows quantum 
trajectories just after the 1st kick. The lines in Fig. 3(a) are trajectories for the single rotor model with no coupling. 
At i = T, 20 particles are settled homogeneously between and 27r on the gi-axis. Figs. 3(b) and 3(c) show typical 
quantum trajectories for the pp-coupling model. For Figs. 3(b) and 3(c), we set 20 particles evenly on the surfaces 
52 = and q2 — 0, respectively, at t = T. The strongly repulsive nature of the trajectory is observed in Fig. 3(a). 
Although the repulsion is also found in Fig. 3(b), for Fig. 3(c) it is partially milder than ones found in Figs. 3(a) and 
3(b). For the coupled rotors, we also see the evolution of the quantum trajectory in the 2-dimensional configuration 
space ((?i,(?2) after the 1st kick. Fig. 4 shows the evolution of the 500 particles which are equally spaced on the 
surfaces q2 = and (72 = at t = T as well as Figs. 3(b) and 3(c). The initial straight row of particles is going to be 
broken at the region where the particle density is increased by condensation of the quantum trajectories. 

Next, we estimate effective acceleration of a quantum particle induced by a kick. In the classical system, nonlinearity 
in the kick is enhanced by the free motion after the kick. For a kicked rotor, therefore, the efficiency of the nonlinearity 
is parameterized hy K = fcT, rather than bare kick strength k. In the quantum system, however, the free motion is 
suppressed by the quantum potential, as seen in Fig. 3. To quantify the effect, we propose to estimate the averaged 
velocity Vi{n) by the distance over which the quantum particle moves during the time interval [nT, (n + 1)T]: 

ViilT-) = ^fe,n+l - 9^n)- (3-3) 

Further, the change rate of the effective velocity Vi is considered as the effective strength of the impulsive kick. 
Therefore, we estimate the effective kick strength by 

= - V,{n - 1)) = ^(((?i,„+i - qi,n) ~ (?»,« - gj,n-i))- (3.4) 



For the classical standard map ( p.2[ ), Fi is approximately a sine function with respect to qi^n- For the quantum 
system, we calculate Fi{n) numerically. At the time t = 30T, when the difference in the dispersion < (pi— < pi >)^ > 
between the single rotor and the coupled ones is apparently observed(Fig.l), we let 500 x 20 probe particles distributed 
uniformly on the qi-coordinate for the single kick rotor, and on the configuration space (91,(72) for the coupled kick 
rotor model. Figs. 5(a) and 5(b) are the results for the single kicked rotor and pp-coupling models, respectively. For 
some particles, the numerical estimation of Vi(29) and Vi(30) does not converge with respect to the integration step 
size, since they may go through the irregularly large quantum potential. In Fig. 5(b), we plot Fi(30) for the 9899 
particles whose deviations in Vi(29) and Vi(30) are smaller than 0.1 if we make the integration step size half. For the 
coupled kick rotors, we can find the tendency that the particles in the region < gi < tt are accelerated and ones in 
the region tt < gi < 27r are decelerated, although not a few particles don't follow the tendency and the acceleration 
and deceleration amplitude is weaker than in the genuine classical case 2.0 sin(gi). The tendency is obscure in Fig. 5(a) 
for a single rotor model. 

We see the evolution of volume elements in the 2-dimensional configuration space (gi, (72)- At t = T, we distribute 
initial probe particles uniformly between and 2n on the surfaces q2 = 2tt x (i=0, 1, • • •, 9) as shown in Fig. 
6(a). Figs. 6(b) and 6(c) show snapshots on the configuration spac e (gi , (72) at t = 2T and t = 3T, respectively. The 
density of quantum particles is conserved in accordance with Eq. ( p. 10 ). The volume elements of the density in the 
configuration space obtain entanglement after several kicks. 

In addition to the dynamics, we also check the distribution of the momentum of the quantum particle, 

fQ{pi) = J j dqidq2\H<li.<l2)?5{p,-^). (3.5) 

By definition, it may differ from the amplitude of the wave function in the pi-representation (Fig. 2). The result is 
given in Fig. 7. 
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IV. SUMMARY AND DISCUSSION 



In order to recognize the role of nonlinearity in quantum systems, we estimated quantum trajectories for the coupled 
kick rotor model. It has recently been known phenomenologically that effect of the nonlinearity appears differently, 
depending on the degrees of freedom of a system [Q- In a single kicked rotor, for example, suppression of the 

momentum diffusion appears and the wave function is localized. For the coupled kick rotors, on the other hand, the 
diffusion in the momentum space seems to continue. In this paper, we compared the single rotor with the coupled 
rotor model from the viewpoint of quantum trajectories. For a single kick rotor, the acceleration of a quantum 
particle induced by the impulsive kicks is suppressed by the quantum potential(Fig. 3(a)). For the coupled two 
rotors, by contrast, entanglement in the two dimensional configuration space is realized (Fig. 6) and suppression of 
the acceleration by the kick is mild in spite of repulsion between trajectories(Fig. 3(c)). In fact, some particles are 
really accelerated (Fig. 5). The dispersion of the momentum of the particles is expanded by the kicks and localization 
is broken for the wave function (Fig. 1). 

A bundle of quantum trajectories flow according to the continuity equation and each trajectory does never go across 
the others due to the uniqueness theory for solutions of differential equations . In one dimensional system, this 
property appears as the phenomenon that the quantum trajectories are repulsive and the complicated density distri- 
bution tends to be flattened. Therefore, the relative acceleration between quantum particles induced by nonlinearity 
is suppressed in the system. This results in the localization of the wave function. For the system with redundant 
degrees of freedom like our coupled kick rotor model, the complicated density distribution is flattened also into the 
second degree of freedom. That is, quantum particles can be scattered forwardly with each other. As a result, the 
repulsive nature due to quantum interference is not straightforward. It is the reason why introduction of the coupling 
yields the possibility that some particles can be accelerated and the localization of the wave function is broken. 

In the de Broglie-Bohm interpretation, the momentum of a quantum particle with a rigid trajectory is equalized to 
the spatial gradient of the phase factor in the wave function. In the pp-coupling model, 

92 = P2 + Cpppi = -q^Vq - Cpp^Vq 

~ 2 dq^R^ dql dql + Oq^dq^ 2 dq^ dqf + Oqj ' + ^""pp'' 

(4.1) 



Equation (4.1) shows that, even if the degrees of freedom, qi and q2 are uncorrelated initially, that is, R{qi,q2) = 
R{qi)R{q2), the coupling induces the gradient of quantum potential with respect to the gi-direction and accelerate 
particles in the direction of the second degree of freedom, q2. This effect decreases the gradient of the density of 
quantum particles and reduces the repulsion in the gi-direction. In the work by Adachi et al 0, the instantaneous 
coupling between the configuration variables qi and q2 is introduced as a model to analyze. In this case, mixing of the 
phase factor in the wave function is induced directly by the coupling and causes the acceleration along the g'2-axis. 

Redundant degrees of freedom help the classical diffusion, which is analogous to the ergodicity although the phase 
space is not compact in our model. Our analysis also shows that multi-dimensional systems emulate partly the 
classical dynamics. On one hand, the quantum trajectory is altered by repulsion in short range. In global range, on 
the other hand, it may express nonlinearity in the system. We may say that the Ehrenfest time scale is not active in 
our analysis where the wave function is spread widely from the beginning and continues diffusion in the uncompact 
phase space. Therefore, there is a possibility that we can treat the appearance of the classical property as a matter 
of the scale of the discrepancy between trajectories rather than a matter of the time scale, such as the Ehrenfest time 
scale. To clarify the dependence of the dynamics on the spatial discrepancy between quantum trajectories, we have 
to make further qualitative analysis. The analysis need more computational memory and more CPU time and will be 
presented elsewhere. 

As demonstrated in this paper, the de Broglie-Bohm picture which gives rigid trajectories is helpful to understand the 
dynamics equivalent to the unitary evolution driven by the Schrodinger equation. In the quantum trajectory picture, 
however, the superposition of the right going and left going waves which have the same and opposite momentum is 
identified with the static state with zero momentum. In our analysis, therefore, we utilized the averaged velocity Vi 



rather than the velocity at a moment, (2.15). In Fig. 2, the characteristic plateaus appears at the classically forbidden 
point {q,p) — (tt, 2n7r)(n — 0, ±1, ±2, ...) which corresponds to the elliptical fixed point, while does not in Fig. 7. For 
comparison with the statistical distribution on the phase space, the Wigner and Husimi descriptions have certainly 
preferable points |p^ . To understand the quantum dynamics in multi-dimensional non-linear systems further, the de 
Broglie-Bohm approach and complementally the Wigner and Husimi functions will be needed. 
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FIG. 1. The time evolution of the dispersion of the momentum, the second order moment Q =< (pi— < pi >)^ > is shown. 
The lines (No) and (P-P) indicate the quantum single kick rotor with the kick strength ki = 2.0 and the kick rotors with the 
momentum coupling, respectively. For the coupled kicked rotors, the kick strength ki and k2 are set to 2.0 and 0.9, respectively, 
and the coupling strength Cpp is 0.2. The initial state is set to the momentum eigenstate with piQ — P2q = n/2 in Eq. (3.1). 



FIG. 2. The momentum distribution | < > [^ at J = lOOT + are shown. The abbreviations (No) and (P-P) indicate 
the single kick rotor and pp-coupling models, respectively. The model parameters ki, k2, and Cpp are the same as in Fig. 1. 



FIG. 3. The quantum trajectory during the time interval from t = T to t = 2T is shown. The horizontal axis indicates the 
time t, and the vertical one shows the gi-coordinate. Fig. 3(a) shows the 20 typical trajectories which are set uniformly along 
the gi-coordinate at t = T for a single rotor. By contrast. Figs. 3(b) and 3(c) show the 20 quantum trajectories which are 
equally spaced on the surfaces 52 = ti" and 52 = 0, respectively, at t = T. The model parameters fci, k2, and Cpp are the same 
as in Fig. 1. 



FIG. 4. For the coupled kicked rotors, the time evolution of probe particles is shown in the time interval from t = T to 
t = 2T. At t — T, 1000 quantum particles are ranged homogeneously on the surfaces g2 = tt and 52 = 0. The horizontal 
and vertical axes show the coordinates gi and 52 of the probes, respectively. Each row like a line shows the snap shot at 
t = LOT, LIT, 1.2T, • ■ ■ , 2.0T. The model parameters fci, k2, and Cpp are the same as in Fig. 1. 



FIG. 5. The effective acceleration _Fi(30) for 500 x 20 probes which are settled a,t t = 30T is shown. Figs. 5(a) and 5(b) 
correspond to the single kick rotor and pp-coupling models, respectively. For the single kick rotor, the probes are distributed 
uniformly between and 27r along the qi-axis. For the coupled kicked rotors, 20 rows which consists of 500 probes along the qi 
axis are arranged uniformly between and 2n. For each probe, the averaged velocities Vi(29) and V2(30) are estimated by the 
two computational runs which are different in the time step size {At = 2.5 X lO-^T and 5.0 x 10"^T). For the coupled kicked 
rotor, the 101 probes whose deviations in Vi(29) and V2(30) are greater than 0.1[1/T] are rejected in Fig. 5(b). The model 
parameters fci, k2, and Cpp are the same as in Fig. 1. 
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FIG. 6. The positions of 5000 x 10 probes axe shown to visualize the deformation of the volume element in the 2-dimen- 
sional configuration space {qi,q2)- Fig. 6(a) shows the initial positions of the probes distributed uniformly on the surface 
q2 = 2-K X i/10(i = 0, 1, • ■ • 9) at i = T. Fig. 6(b) and (c) arc the snap shots for the positions of the probes at t = 2T and 3T, 
respectively. The model parameters fci, fc2, and Cpp are the same as in Fig. 1. 



FIG. 7. The distribution of the momentum of the quantum particles are shown. The distribution function is estimated 
by /q(pi) ~ J J |*&pi5(pi — ^)dgidg2, where $ is the wave function at t = lOOT for the pp-coupling model. The model 
parameters fci, feg, and Cpp are the same as in Fig. 1. 
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